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We describe a new method that is both physically explicable and quantitatively accurate in 
describing the multifractal characteristics of intermittent events based on groupings of rank-ordered 
fluctuations. The generic nature of such rank-ordered spectrum leads it to a natural connection 
with the concept of one-parameter scaling for monofractals. We demonstrate this technique using 
results obtained from a 2D MHD simulation. The calculated spectrum suggests a crossover from 
the near Gaussian characteristics of small amplitude fluctuations to the extreme intermittent state 
of large rare events. 

PACS numbers: 95.30.Qd, 52.35.Ra, 47.27.-i, 94.05. Lk 



Intermittent fluctuations are popularly analyzed using 
the structure function and/or partition function meth- 
ods. These methods investigate the multifractal charac- 
teristics of intermittency based on the statistics of the full 
set of fluctuations. Since most of the observed or simu- 
lated intermittent fluctuations are dominated by fluctu- 
ations with small amplitudes, the subdominant fractal 
characteristics of the minority fluctuations - generally of 
larger amplitudes - are easily masked by those charac- 
terized by the dominant population. It therefore appears 
prudent to search for a procedure that explores the sin- 
gular nature of the subdominant fluctuations by first ap- 
propriately isolating out the minority populations and 
then perform statistical investigations for each of the iso- 
lated populations. Using an example, we demonstrate 
how this idea may be achieved with a rank-order method 
that subdivides the fluctuations into groupings based on 
the range of the scaled-sizes of the fluctuations. 

To be more specific, we consider a generic fluctuating 
temporal event X(t) and form from which a scale depen- 
dent difference series 5X(t, r) = X(t + r) — X(t) for a 
time lag r. We now consider the probability distribution 
functions (PDFs) P(SX,t) of SX(t,r) for different time 
lag values r. If the phenomenon represented by the fluc- 
tuating event X(t) is monofractal i.e., self-similar, the 
PDFs would scale (collapse) onto one scaling function P s 
as follows: 



P(SX,t)t s = P s (6X/t s ) 



(1) 



where s is the scaling exponent. Such one-parameter 
scaling has been suggested for the stock market indices 
[1], magnetic fluctuations [2, 3] and fluctuating events 
of other natural or experimental systems [4]. If the 
PDFs are Gaussian distributions for all time lags similar 
to those characterizing self-similar random diffusion, the 
scaling exponent s is equal to 0.5. For other monofractal 
distributions, the scaling exponent may take on any real 
value. 

In practical applications, expression (1) is sometimes 



approximately satisfied for the full range of the scaling 
variable Y = SX/t s and sometimes only for a portion 
of the range of Y. If the scaling exponent s is obtained 
by estimating the values of the PDFs at SX — 0, the 
PDFs would generally scale at least for a range of Y 
close to the origin. When scaling of the PDFs based 
on (1) is not fully satisfied or only approximately satis- 
fied, the fluctuating phenomenon represented by X(t) is 
multifractal. One conventional method of evaluating the 
degree of multifractal (intermittent) nature of X(t) is to 
study the scaling behavior of the moments of the PDFs 
(conventionally called the structure functions). 



5 m (r) = (\SX(r)\ m ) = 



JX m „ 

\SX(r)\ m P{6X,t)6SX 

where < ... > represents the ensemble average and 5X max 
is the largest value of SX obtainable from the time series 
X{t) for the time lag r. The choice of taking the ensem- 
ble average of the absolute values of the coarse-grained 
differences instead of the values of the raw differences 
is for the purpose of better statistical convergence [5, 
6] . One then proceeds to search for the scaling behavior 
S m (t) ~ . If such scaling is verified for a monofractal 
fluctuating event X(t), the structure function exponents 
would vary linearly with the moment order as £m = Ci m 
where £i = s. If the structure function exponents deviate 
from the above linear relationship, the fluctuating event 
is multifractal. There are several disadvantages of this 
approach. Firstly, the statistical analysis as prescribed 
above incorporates the full set of fluctuations represented 
by X(t). As with most observed PDFs, the statistics are 
generally dominated by those fluctuations with small am- 
plitudes. Thus, the fractal (multifractal) nature of the 
subdominant (larger amplitude) fluctuations is usually 
masked by the fractal nature of the dominant (smaller 
amplitude) fluctuations. Secondly, although deviations 
of the structure function exponents from the linear form 
would indicate that the fluctuating event X(t) is mul- 
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tifractal, the physical interpretation of the multifractal 
nature is not easily deciphered by merely examining the 
curvature of the deviation from linearity. Thirdly, the 
structure functions are usually ill-defined for negative val- 
ues of m. We therefore search for a procedure that would 
remedy the above defects as shown below. 

From the above argument, it appears prudent to per- 
form statistical analyses individually for subsets of the 
fluctuations that characterize the various fractal behav- 
iors within the full multifractal set. We recognize that 
such grouping of fluctuations must depend somehow on 
the sizes of the fluctuations. However, we also realize 
that groupings cannot depend merely on the raw values 
of the sizes of the fluctuations because the ranges will 
be different for different scales (time lags r). Thus, we 
proceed to rank-order the sizes of the fluctuations based 
on the ranges of the scaled variable Y, defined above. 
For each chosen range of AY we shall assume that the 
fluctuations of all time lags will exhibit monofractal be- 
havior and be characterized by a scaling exponent s. The 
question is then how can this procedure be accomplished. 
We continue by constructing the structure functions for 
the chosen grouping of the fluctuations by performing 
moment - structure function - calculations as prescribed 
by (2) with the limits of integral replaced by the end 
points of the range of the chosen AY for each time lag 
r. We then search for the scaling property of the struc- 
ture functions that varies as sm. If such scaling property 
exists, then we have found one region of the multifractal 
spectrum of the fluctuations such that the PDFs in the 
range of AY" collapses onto one scaled PDF. Continuing 
this procedure for all ranges of AY" produces the rank- 
ordered multifractal spectrum s(Y) that we are looking 
for. The determined values of s for each grouping should 
be un-affected by the statistics of other subsets of fluc- 
tuations that are not within the chosen range AY and 
therefore should be quantitatively quite accurate. The 
physical meaning of this spectrum is that the PDFs for 
all time lags collapse onto one master multifractal scaled 
PDF. The spectrum is implicit since Y is a defined as a 
function of s. 

We demonstrate the above outlined procedure via an 
example. The example is based on the results of a large- 
scale 2D magnetohydrodynamic (MHD) simulation. In 
the simulation, ideal compressible MHD equations ex- 
pressed in conservative forms are solved numerically with 
1 024 x 1024 grid points in a doubly periodic (x, z) domain 
of length 2ir in both directions using the WENO code [7] 
so that the total mass, energy, magnetic fluxes and mo- 
menta are conserved. The initial condition consists of 
random magnetic field and velocity with a constant total 
pressure for a high beta plasma. After sufficient elapsed 
time, the system evolves into a set of randomly inter- 
acting multiscale coherent structures exhibiting classical 
aspects of intermittent fluctuations. More detailed de- 
scription of the simulation was reported in one of our 
previous publications [3]. Spatial values (approximately 
one million data points) of the square of the strength of 
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FIG. 1: Raw PDFs as functions of \5B 2 \ for 8 = 8 (solid 
curve in blue), 16 (dotted curve in red), and 64 (dash curve in 
black). |<5-B 2 | is in units of bin size chosen as |5_B 2 |(max)/800. 
There are 800 bins and only first 400 are shown. 
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FIG. 2: Scaled PDFs with s = 0.335. Same line styles and 
colors as in Fig. 1. Y = \8B 2 \/8 S \ \8B 2 \ in units of bin size 
and 8 in grid spacing. 



the magnetic field B 2 are collected over the entire (x, z) 
- plane at a given time and PDFs P(\5B 2 \, S) are con- 
structed for the absolute values of the spatial fluctua- 
tions |<5i? 2 | at scale intervals 6=8, 16, and 64 grid points, 
Fig. 1. Thus, in lieu of temporal fluctuations, the ex- 
ample considers spatial fluctuations. The PDFs are non- 
Gaussian and become more and more heavy-detailed at 
smaller and smaller scales. 

An attempt to collapse the unsealed PDFs according 
to the monofractal scaling formula that is analogous to 
(1) indicate approximate scaling with an estimated scal- 
ing exponent s=0.335, Fig. 2. Structure function cal- 
culations based on the full set of simulated fluctuations 
showed a nonlinear relation between the exponents and 
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FIG. 3: Rank-ordered multifractal spectrum, s(Y). Y = 
\SB 2 \/S S ; \8B 2 \ in units of bin size and S in grid spacing. 
The spectrum is calculated for 8 contiguous ranges of AF. 

the moment order [8]. Because for this example the 
PDFs exhibited approximate monofractal scaling, such 
structure function calculations - though indicating mul- 
tifractality - is strongly masked by the population of the 
smaller fluctuations. Thus, this example represents an 
ideal candidate to test the utility of the new method de- 
scribed in this paper. 

We now proceed to construct the rank-ordered mul- 
tifractal spectrum based on the afore-mentioned proce- 
dure. Thus we sort the fluctuations into ranges of Ay 
between (Yy, F 2 ) with Y = \SB 2 \/S S and evaluate the 
rank-ordered structure functions within each range: 

S m (\SB 2 \,5) = r ' \5B 2 \ m P(\6B 2 \,8)d(\5B 2 \) (3) 

J ai 

where a\ — Y\5 S and a 2 = Y 2 5 s . Expression (3) is a non- 
linear function of s for each moment order to. We now 
search for the value(s) of s such that S m ~ 3 srn within 
each range of the fluctuations so that the rank-ordered 
fluctuations would exhibit monofractal behavior. Inter- 
estingly there exists one and only one value of s in each 
range of AY, that satisfies the above constraint, indicat- 
ing the appropriateness of the ansatz. Unlike the struc- 
ture functions defined for the full range of fluctuations, 
the range-limited structure functions based on (3) exists 
also for negative real values of to. Figure 3 displays the 
calculated rank-ordered spectrum s(Y) based on eight 
contiguous ranges of AY. It is noted that the spectrum 
has values of s ranging between 0.5 and 0.0. The spec- 
trum can be refined by choosing more range intervals 
with smaller range sizes of Ay, although in practice this 
procedure is limited by the availability of simulated data 
points. At Y = 0, the scaling exponent appears to ap- 
proach the self-similar Gaussian value of 0.5. As the value 
of the scaled fluctuation size Y increases the scaling ex- 
ponent decreases accordingly indicating the fluctuations 
are becoming more and more intermittent. At the ex- 
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FIG. 4: Scaled PDFs for Y in (40, 50) for s = 0.15, 0.2, 0.25, 
and 0.3. S —16 (red open circle), 24(black), 32 (blue open 
triangle). Note the triangles are higher than the open circles 
for both s= 0.15 and s=0.25. When s=0.2, they are about 
the same. 

tremely intermittent state, the value of the scaling ex- 
ponent would asymptotically approach the value of zero. 
This would occur at the limit of largest and rarest scaled 
fluctuations. 

For each range of AF, the PDFs would collapse ac- 
cording to its correspondingly calculated exponent value 
s. For example, for the range of Y between (40, 50), the 
PDFs should collapse for the calculated value s « 0.2. 
This is essentially verified as shown in Fig. 4. Figure 
5 shows the results of the rank-ordered structure func- 
tions for the same range of AF. It also indicates scaling 
exponent s = 0.2. 

Such an implicit multifractal spectrum has several ad- 
vantages over the results obtainable using the conven- 
tional structure function and/or partition function cal- 
culations. Firstly, the utility of the spectrum is to fully 
collapse the unsealed PDFs. Secondly, the physical in- 
terpretation is clear. It indicates how intermittent (in 
terms of the value of s) are the scaled fluctuations once 
the value of F is given. Thirdly, the determination of the 
values of the fractal nature of the grouped fluctuations is 
not affected by the statistics of other fluctuations that do 
not exhibit the same fractal characteristics. Fourthly, it 
provides a natural connection between the one-parameter 
scaling idea (1) and the multifractal behavior of intermit- 
tency. 

To summarize, we have introduced a new rank-ordered 
procedure based on the sizes of the scaled fluctuations 
to view the multifractal nature of intermittent fluctua- 
tions. The suggested implicit multifractal spectrum anal- 
ysis provides a physically meaningful description of inter- 
mittency and is quantitatively accurate because of the 
cleanliness of the procedure of statistical sampling. The 
method can easily be generalized to situations of higher 
dimensions as well as correlation and response functions 
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of several independent variables involving intermittency 
of spatiotcmporal fluctuations. 
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FIG. 5: Results of rank-ordered structure functions for Y in 
(40, 50) for s = 0.15 (circle), 0.2 (star), and 0.25 (triangle). 
Cm is defined by S m (\SB 2 |, 5) ~ 6^ m and is obtained by a 
linear fit of \og(S m (\5B 2 \, 5)) versus log((5) for 3 from 16 to 32. 
The dash line, solid line and dotted line show linear functions 
srn for s=0.15, 0.2, and 0.25, respectively. For s=0.2, £ m « 
sin. 
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